Photo-induced states in a Mott insulator 
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We investigate the properties of the metallic state obtained by photo-doping carriers into a Mott 
insulator. In a strongly interacting system, these carriers have a long life-time, so that they can 
dissipate their kinetic energy to a phonon bath. In the relaxed state, the scattering rate saturates at 
a non-zero temperature-independent value, and the momentum-resolved spectral function features 
broad bands which differ from the well-defined quasi-particle bands of a chemically doped system. 
Our results indicate that a photo-doped Mott insulator behaves as a bad metal, in which strong 
scattering between doublons and holes inhibits Fermi-liquid behavior down to low temperature. 



Carrier doping provides a convenient way of control- 
ling the properties of strongly correlated materials. The 
high-T c cuprates, e.g., undergo a doping-driven transi- 
tion from an antifcrromagnetic Mott insulator to a cor- 
related metal that becomes superconducting at low tem- 
peratures. While doping is usually achieved by modify- 
ing the chemical composition of a material, photo-doping, 
i.e., a change of the electron and hole concentrations by 
irradiating light on a sample, provides a straightforward 
way to influence material properties in an ultra-fast man- 
ner. Photo-induced insulator-to-metal transitions have 
been demonstrated in organic materials [l], charge- 
transfer insulators 0] and cuprates 

In contrast to chemically doped states, photo-doped 
states simultaneously have electron and hole-like carriers, 
and these carriers are typically inserted with large kinetic 
energy. In a semiconductor, this mainly results in a dif- 
ferent occupation of valence and conduction bands, which 
themselves are rigid, i.e., independent of filling and tem- 
perature. In contrast, the formation of a Fermi liquid in 
a doped correlated system leads to a narrow band of de- 
localized quasi-particle states without any counterpart in 
the insulator. Because these quasi-particles are strongly 
damped at high temperature, and because the adiabatic 
correspondence between bare electrons and quasiparticle 
excitations of a Fermi liquid is anyway difficult to rec- 
oncile with the presence of both holes and electrons, the 
properties of a correlated photo-doped state are not ob- 
vious. Its large kinetic energy works against correlations 
and could cause a rigid-band behavior. On the other 
hand, life-times of a few picoseconds indicate that carri- 
ers are strongly coupled to spin fluctuations or phonons 
d, 3 . They might thus dissipate their energy before re- 
combination, and reveal correlation effects in a photo- 
doped state and possible differences to a Fermi liquid. 

In this paper we contrast photo-doped and chemically 
doped Mott insulators in the simplest possible setup, 
which is a paramagnetic single-band Mott insulator that 
is initially perturbed by an intense laser pulse. For this 



purpose we focus on the Hubbard model, 

which describes electrons that hop between nearest 
neighbors (ij) of a crystal lattice and interact through a 
local Coulomb repulsion U. The electric field of the laser 
pump is determined by the vector potential, E = —^d t A, 
which in turn enters Eq. ([T]) by the Peierls substitution, 
i.e., band energies are shifted to Ik = e k-ea/hcA(t) ( see > 
e.g., Ref. Q). If the Mott gap is large (U ^ v), electron- 
hole pairs have a long lifetime, because the electronic 
decay-channel involving a direct transfer of interaction 
energy into kinetic energy becomes inefficient Q . To ac- 
count for dissipation of energy to other degrees of free- 
dom, we weakly couple our model to a bosonic bath. As 
we will show below, the resulting reduction of the kinetic 
energy of the particles drives the system into a metallic 
state that is not accessible in thermal equilibrium. 

The model is treated within the noncquilibrium ex- 
tension 0, HI of dynamical mean-field theory (DMFT) 
A direct way to incorporate dissipation into this 
formalism is to attach one particle-reservoir to each lat- 
tice site [IB-diJ . Here we instead use baths of harmonic 
oscillators, such that the particle number remains con- 
stant while energy is exchanged with the environment. 
The electronic self-energy, which is a functional of the 
local Green function G in DMFT, is then the sum of 
the electronic contribution Eu[G], and a bath contribu- 
tion. Vertex corrections are neglected for weak coupling 
to the bath. The bath contribution is the lowest order 
diagram for a Holstein-type electron-phonon coupling, 
Ediss[G] = X 2 G(t,t')D(t,t'), where A measures the cou- 
pling strength, and D(t, t') is the equilibrium propagator 
for a boson with energy uq at given temperature 1/(3 
(such that the bath has no memory); following the nota- 
tion of Ref. [lH for contour-ordered Keldysh Green func- 
tions, D(t,t') = -iTr[T c exp(-i J c dtw tfb)b(i)tf(t')]/Z 
(luq = v in the following). 

Local diagrams for Ey [G] are summed to all orders by 
solving an auxiliary single-impurity Anderson problem, 
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just as it is done for DMFT without dissipation. Below 
we will use the self-consistent strong-coupling expansion 
[T3 | to solve this impurity model. Because dissipation 
introduces a timescale oc 1/A 2 to the dynamics, which 
may be much longer than the hopping time (and which 
we must access within our simulation), we are restricted 
to the lowest order of the strong-coupling expansion, the 
non-crossing approximation (NCA) [To] ]. However, our 
main results are obtained for a (doped) Mott insulator 
at U 3> v, above and close to the Fermi-liquid coherence 
temperature, where the use of NCA is still justified (U = 
Uv, as for the organic salt ET-F 2 TCNQh The resulting 
equations have been explained in Ref. [16J. To include 
dissipation as described above we just have to replace 
the band energy l\~ by ik + Sdiss[G] in those equations. 

DMFT is exact for infinite dimensions [l7| and yields 
the generic behavior of a high-dimensional lattice, in- 
dependent of the lattice geometry used within the cal- 
culation (mainly, the energy scales are renormalized for 
different geometries) . In this paper we perform all calcu- 
lations for the paramagnetic phase of a one-dimensional 
system, because this drastically simplifies the momen- 
tum summations, but we expect the results to be repre- 
sentative of high-dimensional systems. The unit of en- 
ergy is set by the hopping v, and time is measured in 
units of h/v. In these units, the critical end-point of 
the first-order Mott transition line for the half-filled sys- 
tem is located at U « 4-5. The pump pulse, E(t) = 
E sm[Q(t — 1 )] exp[-4.605 X (i — £o) 2 /*o]j i s taken to be 
a sine- wave with frequency fi, amplitude Eq, and Gaus- 
sian envelope. The number of cycles, N = to^lf-K, is 
taken between 2 and 10 below. To analyze photo-doped 
and equilibrium states, we compute the real-time opti- 
cal conductivity [l8| and its partial Fourier transform, 
a(u!,t) = f Q ma *ds cr(t,t — s)e lujs . We also define an aver- 
age of <r{u!,t) which is more local in time and frequency 
(and almost independent of the cutoff s max = t), 

a(w, t; A) = f "^ds e~ s2A ^ 2 a{t, t - s)e ius , (2) 
Jo 

equivalent to broadening a(ui,t) with exp(— w 2 /2A 2 ). 

Before studying the photo-doped case, we investigate 
the influence of the dissipative bath on equilibrium prop- 
erties of the (doped) Mott-insulator. Pronounced bath- 
induced artifacts begin to appear for A > 0.8, as exem- 
plified for the spectral function of the insulating phase 
in Fig. QJ,. Because our approximation is only suitable 
for weak coupling between electrons and bath, we will 
from now on restrict A to values where those features 
are absent and dissipation implies only a slight broad- 
ening of the spectrum (A < 0.6). Small electron doping 
(n > 1) of the Mott insulator leads to the formation of a 
narrow quasi-particle peak at the lower edge of the up- 
per Hubbard band, which goes along with the emergence 
of a Drude peak in the low-frequency part of the opti- 
cal conductivity (Fig. [Tb), a{w) = D^ir~ 1 [y 2 + lo 2 }^ 1 + 
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FIG. 1: (a) Equilibrium spectral function of the half-filled 
insulator for U = 8, /3 = 10, and various A. (b) Frequency- 
dependent conductivity cr(oj) of the doped Mott insulator at 
various temperatures, U = 14, n = 1.02. Fourier artifacts due 
to a finite cutoff s max = 40 appear when the peak becomes too 
narrow (/? = 15). The inset shows the Drude relaxation time 
I/7, obtained from a Lorentzian fit to o(uj) for < uj < 0.2. 

Cincoh(w). For a Fermi liquid, the scattering rate 7 de- 
creases with decreasing temperature, with the asymp- 
totic behavior 7 ~ T 2 . Although the NCA approxima- 
tion is known to suffer from non-causal artifacts deep in 
the Fermi liquid regime, the onset of this Fermi liquid 
regime can still be observed around (3 = 20, by extract- 
ing 7 from a fit of ct(qj) with a Lorentz curve (Fig. [TJ 
inset). For /3 < 20, coupling to a bath with A = 0.6 has 
only a small influence on the scattering rate. 

In the following paragraphs we study the relaxation 
of the double occupancy and kinetic energy (Fig. [2ji-[2f ) 
after a pump at f2 = U has created a small number of 
charge excitations in the system. Without dissipation 
(A = 0) and for small U (U = 4, 5), both d(t) and E kin (t) 
follow an exponential relaxation d(t) ~ doo + Ae~'/ T to- 
wards a value d^ ^ d(0). This fact has already been 
described for the hyper-cubic lattice in Ref. where 
it was also verified that doo is consistent with a thermal- 
ization of the system at a temperature that is determined 
by the total amount of absorbed energy. The thermal- 
ization time exponentially increases with U (Fig. [2jr), in 
agreement with the predicted exponentially long lifetime 
of doublons in the Hubbard model for U > w In 
the insulator (U = 8), r becomes so large that d(t) and 
E]u-a(t) remain constant on the scale of the plot, apart 
from a short initial transient. 

Dissipation completely modifies this picture: Instead 
of approaching a thermal equilibrium value at higher 
temperature, both d(t) and -Ekin(i) now relax back to 
their initial values for U = 4 and U = 5. For [7 = 4, this 
process can be understood within a two-temperature pic- 
ture [l9[ : Because the electronic thermalization time is 
only a few h/v (Fig. 0i), electrons can reach a quasi- 
equilibrium state at high temperature T* before energy 
is transferred to the environment. The temperature T* 
is subsequently reduced towards the fixed temperature of 
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FIG. 2: Relaxation of kinetic energy (a-c) and double occu- 
pancy (d-f) for U = 4 (left), (7 = 5 (middle), and U = 8 
(right), during and after a pump with frequency Q = U, 
N = 10 cycles, i?o = 1; /3 = 10. Curves show the difference 
between the time-evolved values and the initial state values, 
brought to the same scale by an additional factor ("x"). Bold 
solid lines correspond to no dissipation, (g) Same data as (d), 
plotted against rescaled time t' = X 2 (t — to), where td is the 
duration of the pump. The inset shows the equilibrium value 
of d as a function of /3 for U = 4. (h) Relaxation time of the 
double occupancy, from an exponential fit. 



the bath with a rate entirely determined by A. This inter- 
pretation is supported by the observation that the curves 
d(t) — d(0) for [7 = 4 fall on top of each other when the 
time axis is rescaled by A 2 (Fig. [2g). indicating a pas- 
sage through the same sequence of equilibrium states. 
In fact, because U = 4 is in the metal-insulator crossover 
region, the double occupancy behaves non-monotonically 
as a function of temperature (inset of Fig. [2g) , which ex- 
plains the minimum in d(t). 

Because the energy quantum that can be transferred to 
the bath is limited (w = v), coupling to the bosons does 
not open an efficient channel for the doublon decay for 
U ^> v. Consequently, the relaxation rate of the double 
occupancy to its initial value shows a similar exponen- 
tial increase with U as the thermalization time for A = 
(Fig. (2Ji), and the bath has no effect on d(t) in the insu- 
lator (U = 8). The kinetic energy, however, does relax to 
a value that is lower than after the excitation, indicating 
that doublons and holes are cooled to the temperature of 
the bath before they recombine. In the following we will 
contrast the resulting low-energy state with a chemically 
doped metal, by comparing their conductivities. 

Figure shows the optical conductivity of the photo- 
doped Mott insulator for various time-delays after the 
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FIG. 3: (a) Optical conductivity <r(u),t) [Eq. ©, A = 0.2] 
at various times for the photo-doped Mott insulator (U = 14, 
/3 = 6, A = 0.6; Pump: Ea = 2.5). (b) Scattering time as a 
function of (bath) temperature for the chemically doped case 
(/i = 4.6) and the photo-doped case. Pump: Q. = 14; exci- 
tation density 8 w 0.8% {E = 1.5) and 8 w 2% (E = 2.5). 
Circle symbols: Eq. ((3]) taking a(u},t) at t = 70. Diamonds: 
Eq. (|3} taking a{uj,t) from the extrapolation to t = oo. Ar- 
row: 7a for a temperature quench (8 « 2%), see text. (Inset) 
a(u>,t) — a(aj,oo) for ui — 0, where ct(oj,oo) is obtained from 
a fit a(ui, t) = (j{uj, oo) + Aexp(— t/r). 



excitation. The decrease of the kinetic energy due to dis- 
sipation leads to an increase of spectral weight accord- 
ing to the sum rule Jduj a(uj,t) = —E\an(t)- Nearly all 
added weight enters the Drude peak, while the weight 
W^Hub in the Hubbard band is reduced with respect to 
the undoped case. In fact, W^kub measures only near- 
est neighbor short-time correlations, and with respect to 
this, photo-doped and chemically-doped states just look 
alike: the reduction of Wnub is time-independent and 
proportional to doping S = Ad + Ah, where here and in 
the following Ad and Ah denote the differences in the 
doublon and hole densities with respect to the un-doped 
insulator (6 w 2% in Fig.©. 

The Drude peak and its width 7, on the other hand, 
can give a hint on whether or not coherent quasi-particles 
are being formed. Instead of 7 we compute the approxi- 
mate measure 



7A = uJy / a{uj = 0;A)/a{u>;A)-l, 



(3) 



which can be obtained from a finite time-window (after 
relaxation of -Ekm), and approaches 7 for A — > and a 
Lorentz curve (we choose ui = A = 0.5). Although the 
asymptotic behavior 7 oc T 2 is cut off below the scale A, 
7a still reveals the increase of the inverse scattering rate 
in the interesting temperature range where the chemi- 
cally doped state enters the Fermi liquid regime (Fig.[3]D). 
In the photo-doped state, in contrast, I/7A saturates at a 
smaller value for T < 1/5. Note that a is not yet fully sta- 
tionary at the largest simulation time t = 70. However, 
an exponential extrapolation, a(t) = a(oo)+Acxp(—t/T) 
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(Fig. [3l inset) corrects 7a only slightly, and moreover, it 
leads to a temperature-independent value (diamond sym- 
bols in Fig. (SJd). This finding, which is a central result 
of this paper, is rather insensitive to the excitation den- 
sity (see results for E = 1.5, 8 « 0.8%, and E = 2.5, 
8 w 1 — 2%, in Fig. [5b). It shows that the photo-doped 
state does not behave like a Fermi-liquid upon lowering 
the temperature. 

Further differences between photo-doped and chemi- 
cally doped states are evident from their spectral func- 
tions (Fig. 2}: The occupied density of states A < (u,t) = 
(l/-7r)Re J °° dse luJS (c* (t+s)c(t)) , which for a quasi-steady 
state is related to a time-resolved photocmission spec- 
trum j20l ]. evolves from a broad photo- induced distri- 
bution for small times into a peak concentrated at low 
frequencies. At the same time, the density of states, 
A(u,t) = (l/V)Re dse luJS ({c(t+s)\c(t)}), is not rigid 
as it would be for a noninteracting system, but it devel- 
ops a feature at the lower edge of the upper Hubbard 
band (the lower Hubbard band is symmetric). This fea- 
ture is different from the quasi-particle peak of a chemi- 
cally doped system, both when the density n of the latter 
is adjusted to the total number of photo-doped carriers 
(n — 1 = 2Ad) or to the number of electron-like carriers 
(n — 1 = Ad). A closer look at the momentum resolved 
spectral function A k (uj,t) = (l/7r)Re / °° dse luls ({c k (t + 
s)',Cfc(f)}) after the system has become nearly station- 
ary (t = 40) reveals that this feature belongs to a rela- 
tively broad band of heavily scattered charge excitations 
(Fig. |H inset). 

For U ^> v, strong-coupling perturbation theory for 
the Hubbard model gives a possible interpretation of 
our results: Up to corrections of higher order in v/U, 
a Schricffcr- Wolff transformation brings the Hamiltonian 
([T]) into a form Hu^ v which separately conserves both 
doublon and hole numbers [2lj]. Projecting Hjj^ v onto 
the low-energy subspace without doublons (for the hole- 
doped case) or without holes (for the electron-doped 
case), leads to the familiar t- J- model. In nonequilibrium, 
however, interesting dynamics can occur in different sec- 
tors of the Hilbert space, such as metastable supercon- 
ductivity in the absence of singly-occupied states [Hj], 
or the strong-coupling prethermalization after an inter- 
action quench |23j. The simplest possible photo-doped 
state after relaxation is a low-temperature state of Hjj^ v 
with equal density of both doublons and holes. Our re- 
sults hence suggest that this corresponding equilibrium 
state is a bad metal rather than a Fermi liquid in the 
investigated temperature range, in contrast to the case 
of pure electron or hole doping. At present we cannot di- 
rectly compute equilibrium properties of £/(/>„ for both 
electron and hole doping. However, we can verify that the 
photo-induced state is rather independent of the initial 
preparation, which is a prerequisite for being a thermal 
state: We have performed a similar numerical analysis as 
above, but now the system is initially prepared at high 
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FIG. 4: Spectral function in the upper Hubbard band. 
A(u,t), A < (u),t): photo-doped, U = 14, /3 = 15, A = 0.6; 
pump n = 14, Ea = 2.5 (excitation density Ad = 0.01). 
Curves n = 1.01 = 1 + Ad and n = 1.02 = 1 + 2Ad cor- 
respond to chemically doped states, where frequency axis is 
offset by Aoj = fi. Inset: Momentum-resolved Ak{ui) for var- 
ious k and same parameters. Light red curves: chemically 
doped case at n = 1.02; shaded curves: photo-doped case. 



temperature ((3* = 0.4), and then suddenly coupled to a 
bath at lower temperature (j3 = 10). In this somewhat 
artificial temperature-quench, carriers in the system are 
of thermal origin rather than photo-excited, but the long- 
time behavior of 7a nevertheless turns out to be the same 
as for the photo-doped case (arrow in Fig. \3jp). 

In conclusion, we have demonstrated that the prop- 
erties of photo-doped and chemically doped Mott insu- 
lators are quite distinct, even if photo- induced carriers 
are allowed to transfer their initially high kinetic en- 
ergy to the environment before recombination occurs: A 
temperature-independent scattering time of the carriers 
in the photo-doped system indicates that when a Mott 
insulator is doped with electrons and holes at the same 
time, its properties remain non-Fermi liquid like down to 
much lower temperatures than for pure electron or hole 
doping. Experimentally, this results in a low mobility of 
photo-doped carries, which could be seen in a THz-study 
of the Drude part of photo-doped metals. On a more 
fundamental level, our work raises questions concerning 
the possible ground states of the limiting strong-coupling 
model for the Hubbard model with both electron and hole 
doping. Does Fermi liquid behavior appear below a new, 
strongly reduced, coherence scale? Should we expect new 
quantum phases formed from paired doublons and holes? 
Or will this low-energy state simply phase-separate into 
purely hole-doped and electron-doped regions? 
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